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Abstract 

We review the concept of r-function for simple analytic curves. The r- function 
gives a formal solution to the 2D inverse potential problem and appears as the 
r-function of the integrable hierarchy which describes conformal maps of simply- 
connected domains bounded by analytic curves to the unit disk. The r-function 
also emerges in the context of topological gravity and enjoys an interpretation as 
a large A'^ limit of the normal matrix model. 



1. Recently, it has been realized ||T], H that conformal maps exhibit an integrable struc- 
ture: conformal maps of compact simply connected domains bounded by analytic curves 
provide a solution to the dispersionless limit of the 2D Toda hierarchy. As is well known 
from the theory of solitons, solutions of an integrable hierarchy are represented by r- 
functions. The dispersionless limit of the r-function emerges as a natural object associ- 
ated with the curves. In this paper we discuss the r-function for simple analytic curves 
and its connection to the inverse potential problem, area preserving diffiomorphisms, the 
Dirichlet boundary problem, and matrix models. 



2. Inverse potential problem. Consider a closed analytic curveQ 7 in the complex plane and 
denote by and D_ the interior and exterior domains with respect to the curve. The 
point 2; = is assumed to be in Z)+. Assume that the domain is filled homogeneously 
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with electric charge, with a density which we set to be equal to 1. The potential $ 
created by the charge obeys the equation 



1 if z = X + iy E 
if z = X + iy E D_ 



The potential $ can be written as an integral over the domain D^: 

^{z,z) = -- f <fz' log|z-/| (2) 

TT JD+ 

In the exterior domain D_, the potential is the harmonic function whose asymptotic 
expansion as 2; — > cxd is given by 

$-(^,^) = -2tolog|2| +27^e^^2-^ (3) 

fc>0 ^ 



where 

Vk 



IT 

are the harmonic moments of the interior domain and 



- / z'^dh {k > 0) (4) 

TT Jd+ 



nto = / d^z (5) 



D 

is its area. In the interior domain D^, the potential (Q) is equal to a function which 
is harmonic up to the term — The expansion of this function around z = is 

$+(2,^) = -^;o + 27^e^^fc2^ (6) 

fc>0 

Here 

tk = -\ I z-'dh (k > 0) (7) 

TTK JD- 

are the harmonic moments of the exterior domain and 
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Vo = - 
TT 



/ \og\z\d'^z. 
Jd+ 



The two sets of moments and (0) are related by the conditions that = $ , 
9^$"*" = dz^~ on the curve 7. 

The inverse potential problem is to determine the form of the curve 7 given one 
of the functions $^ or $~, i.e. given one of the infinite sets of moments. We will 
choose as independent variables the area vrto and the moments of the exterior tk {k > 1). 
Under certain conditions, they completely determine the form of the curve as well as the 
moments Vk {k >0) 0. More precisely, {tfcj^o ^ good set of local coordinates in the 
space of analytic curves. For simplicity we assume in this paper that only a finite number 
of tfc are non-zero. In this case the series @ is a polynomial in z, z and, therefore, it 
gives the function for z G -D+. Note that t^, vq are real quantities while all other 
moments are in general complex variables. 
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3. Variational principle. Consider the energy functional describing a charge with a 
density p{z, z) in the background potential created by the homogeneously distributed 
charge with the density +1 inside the domain (|T]): 

£{p} = J J (fz(fz' p{z,z) log \z — z'\p{z' , z) JcPz p{z,z) ^{z,z). (9) 

The first term is the 2D "Coulomb" energy of the charge while the second one is the 
energy due to the background charge. Clearly, the distribution of the charge neutralizing 
the background charge gives the minimum to the functional: po = —1 inside the domain 
and po = outside. At the minimum the functional is equal to minus electrostatic energy 
E of the background charge : 

- E = mmp£{p} = ^ f (fz [ Sz! log \z - z'\ = — - [ dPz $(2, z). (10) 

TT^ JD+ JD+ ZTT JD+ 

Varying over p and then setting p = — 1 inside the domain, we obtain eq.(§). 

The first corollary of the variational principle is that the E is a. potential function for 
the moments. Eq. (||) suggests to treat vq and tk as independent variables, so moments of 
the interior, Vk, k > 1, and to are functions of vq and tk- Let us differentiate E or —S{p} 
at the extremum with respect to the parameters Vq, tk- Since po minimizes the functional, 
the derivative is equivalent to the partial derivative of S at the fixed extremum p. This 
gives 

dE dE _ dE 

dtk dtk ovo 

where the partial derivative with respect to tk is taken at fixed Vq and tj [j 7^ 0,k). 
Therefore the differential dE reads 

dE = ^{vkdtk + Vkdtk) - todvQ. (12) 

fc>0 

Let us note that the variational principle may be formulated in a number of different 
ways. One particular variational principle is suggested by the matrix model discussed in 
the Sec. 9. In this case one consider a charged liquid in the potential 

V{z,z) = zz + Vo-Y.{ikZ^ + tkz'') (13) 

A;>0 

defined everywhere on the plane and vq and tk are parameters. The energy of the charged 
liquid 

£{p, 1/} = _i_ fSzj d^z p{z, z) p{z', z) log\z- z'\ +- f d^z p{z, z) V{z, z). (14) 



reaches the minimum if the liquid forms a drop with the density po = — 1 bounded by 
the curve determined by parameters of the potential vo and tk- For another version of 
the variational principle see 0. 

4. T-function- It is more natural to treat the total charge to rather than vq as an 
independent variable, i.e. to consider the variational principle at a fixed total charge 
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/ log 








Jd+ 


z 





to = J p(Pz. This is achieved via the Legendre transformation. Let us introduce the 
function F = E + toVo , whose differential is 

dF = ^{vkdtk + Vkdtk) + vodto. (15) 

k>0 

We define the r-function as r = e^, so that 

' ' dhdh'. (16) 

The r-function is a real function of the moments {to, ti, ^2, • • •}• Under the as- 
sumption that only a finite number of them are non-zero, we can substitute (|^) into 
(|T6|) and perform the term-wise integration. Taking into account that ^ \z\'^d'^z = 
1^0 + \ ^k>ok{tkVk + ikVk) (a simple consequence of the Stokes formula), we get the 
expression for the r-function in terms of tk and Vk'- 

21ogr = -Ul + tovo -\Y.^k- 2){tkVk + ikVk). (17) 

^ ^ fc>0 

Rephraising (0) we get the main property of the r-function, which has been used as 
its definition in Ref.0 

(91ogr 91ogr _ 51ogr 

— = Vk, ^ = Vk, -jr. — = ^^0 18 
otk dtk dto 

where the derivative with respect to tk is taken at fixed tj {j ^ k). 

Two immediate consequences of the very existence of the potential function are sym- 
metry relations for the moments 

dvk _ dvn dvk_ _ dvn 
dtn dtk ' din dtk 

and the quasi-homogeneity condition for the r-function: 



(19) 



4 log r = -4 + 2(„^ - - 2) {fj-^ + tj-^) ^ (20) 



Apart from the term —t\^ this formula reflects the scaling of moments as z — > \z: 

tk ^ A2-4fc (A; > 0), Vk ^ X^^'^Vk {h>\). 

As an illustration we present the r-function of ellipse 0. In this case only the first 
two moments t\ and ^2 are nonzero[|: 

logr = -^tl + -tn log I ^° I H ^V-TT {\t\? + t\t2 + t?t2 



^The r-fuction for the ellipse (at t\ — 0) appeared in Rcf. |^ as the limit of the Laughlin wave function 
or a planar limit of the free energy of normal matrix models, see Sec. 9 



4 



5. Schwarz function and generating function of the conformal map. Consider a univalent 
conformal map of the exterior domain D_ to the exterior of the unit disk and expand it 
in Laurent series: 

w{z) = -z + Y^PjZ-^, (21) 

j=o 

where the coefficient r is chosen to be real and positive. The series for the inverse map 
(from the exterior of the unit disk to D_) has a similar form: 

oo 

z{w) = rw + ^ UjW^^ . (22) 

j=0 

Chosen w on the unit circle, eg . (12^) gives a parametrization of the curve. By the definition 
of an analytic curve, the map can be analytically continued to a strip-like neighborhood of 
the curve belonging to The continuation is given by the Riemann- Schwarz refiection 
principle (see e.g.0): 

w = iw{S{z)))-\ (23) 

where S{z) is the point refiected relative to the curve^ Following 0], we call S{z) the 
Schwarz function of the curve. Let us recall its construction. Write the equation for the 
curve F{x,y) = in complex coordinates, F(^^, ^^) = 0, and solve it with respect 
to z. One gets the Schwarz function: z = S{z). The Schwarz function is analytic in a 
strip-like domain that includes the curve. On the curve the Schwarz function is equal 
to the complex conjugate argument. The main property of the Schwarz function is the 
obvious but important unitarity condition 

S{S{z)) = z (24) 

(the inverse function coincides with the complex conjugate function). In terms of a 
conformal map the Schwarz function is 

oo 

S{z) = rw'\z) + J2 UjW^{z). (25) 

j=0 

Using the Schwarz function one can write the moments of the exterior and the interior 
domains (||,^ as contour integrals 

tn = ^{ z-''S{z)dz, v^ = ^cf z^S{z)dz (26) 

Eq. (PBD yields the Laurent expansion of the Schwarz function 

oo J oo 

S{z) = J2 ktkz"-' + - + E ^kz-"-'. (27) 

k=l ^ k=l 

Now let us define the generating function Q{z), related to the Schwarz function by 

S{z) = d,Q{z)- (28) 



■^We use the notation: given an analytic function f{z) — J^j fj^''i set f{z) = J^j fj^''- 
''This is due to a more general statement f{z)(Pz = ±^ §^ f{z)S{z)dz, where f{z) is an analytic 
function in the domain D±. 
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The latter is given, according to (0), by the Laurent series 

OO 1 oo 

Q{z) = Y,tuz''--v, + to\ogz-Y.^z-^ (29) 
k=i ^ k=\ ^ 

It can be represented as = ^^^\z) + — ^fo, where are analytic in 

D± respectively: 

n^^\z) = - / log(l - ^)d'z' = J2tkz' (30) 



k=l 

oo 



n^-\z) = - [ \ogiz - z')dh' = tologz - f; ^z-'^ (31) 

From (IJI) we see that ^-{z,z) = -2nen^-\z) and $+(z,^) = 2neQ^+\z) - vq - \z\'^. 
Contrary to the potentials "l*^, the analytical functions VL^ and —VL~ do not match each 
other on the curve. The discontinuity gives the value of the generating function restricted 
to the curve 

VL{z) = -\z\^ + 2iA{z), ze-i (32) 

where A{z) is the area of the interior domain bound by the ray Lp = aigz and the real 
axis. As a corollary, it is easy to show that variations of the Q{z) on the curve with 
respect to the rea/ parameters to, TZetk and Xmtk are purely imaginary. This allows one 
to apply the Riemann-Schwarz reflection principle to analytical continuation of 

Hk{z) = dt,n{z), Hk{z) = ~\Q{^) (33) 

and to prove the fundamental relations 

duM^) = \ogw{z), (34) 

d,M^) = {z''{w))^ + l{z'^{w))^ (35) 

d,M^) = {S'^{z{w)))_ + ^-{S\z{w)))^ (36) 

The symbols {f{w))± mean a truncated Laurent series, where only terms with positive 
(negative) powers of w are kept, while {f{w))o is the constant term (w°) of the series. 
Note that the derivatives in eqs. (|3^|36D are taken at fixed z. 

To prove (p4|), we first notice that 

dto^{z{w)) = \ogz — ^^-^+negative powers in ^ = logwr ^^-^ + negative powers in w. 

Independently, one can show that dt^vo = 21ogr. 

Then, using the Riemann-Schwarz reflection principle, we may write dtQ^l{z{w)) also 
in the form dtgCl{S{z{w)). Expanding the latter in S{z) and then, using expansion of 
(P5|) in w, we have 

dtQ^{S{z{w)) = log5'(z) [-negative powers in S{z) = log w-|-positive powers in w. 
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Comparing both expansions, we conclude that dtf^Q{z) = logw{z). Similar arguments 



are used in the proof of (^) and (36) 



6. Dispersionless Hirota equation and the Dirichlet boundary problem. Using the rep- 
resentation (|18D of the moments Vk as derivatives of the r-function, one can express the 
conformal map w{z) (|^) through the r-function: 

\ogw = \ogz- dto (^dto + -^"^tfe j log^- (37) 
With the help of the r-function, eqs. (|35|j36|) can be similarly encoded as follows: 

dA^og{wiz)-wiC)) = ^^^,+ \T.^''''\] (Er"~X) logr (38) 



Kk>l / \n>l 



- d.d^\og{w{z)w{C) - l) = Er"-'5t„ logr (39) 

\fc>l / \n>l / 

The derivation is similar to the one given in P, |^ for the case of the KP hierarchy. 
Moreover, these equations in the integrated form are most conveniently written in terms 
of the differential operators 

m = E '-^d,, , D{z) = E (40) 
fe>i k>i 



From (38,39) one obtains: 

w{ z) - w{C) 

z-C 2 



log ' i = --dl\ogr + Diz)DiC)\ogT (41) 



- log (l - , ^.-. ] = D{z)D{0\ogT (42) 

Combining (^) and (^l)) , one obtains the dispersionless Hirota equation (or the dis- 
persionless Fay identity) for 2D Toda lattice hierarchy |0: 

Eq. (^3l), after being expanded in powers of z and C, generates an infinite set of relations 
between the second derivatives c^t^^t^logr of the r-function. Using (|^) instead of (|4T|), 
a similar equation for the mixed derivatives dt^di^logr can be written: 

Let us conclude this section with two other forms of the dispersionless Hirota equation 
for the conformal map. They emphasize a relation between the Hirota equation and two 
fundamental objects of the classical analysis: the Green function of the Dirichlet problem^ 
and the Schwarz derivative. 



^This relation is pointed out to us by L. Takhtajan. 
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The Green function of the Dirichlet boundary problem for the Laplace operator in 
D_ expressed through the conformal map w{z) is: 



G(z,C)= log 



w{z) - w{Cj 



w{z)w{() — 1 



(45) 



Combining (^Tj) and (^21), and using the notation P0|), we represent the Green function 
as follows: 

2G{z, = 2 log \z-' - C'\ + (5*0 + D{^) + D{z)) + D{0 + D{0) logr (46) 

This formula generalizes ( pT] ) since (^) becomes the real part of (0) as C oo. (As 
C — i> oo, G{z,Q ~^ ~^og\w{z)\.) Note also that the real part of (|37D can be written in 
the form 

$(2, z) = -2to \og\z\ + (d{z) + D{z)) log r 

where $ is the potential (H) {z G DJ). 

The l.h.s. of eq.(^) generalizes the Schwarz derivative of the conformal map 



T{z) 



w"'{z) _ 3 f w"{zY 
w'{z) 2 \ w'{z) _ 



6 lim^^^(9^9^1og 



w{z) — w{Q 



(47) 



Taking the limit ^ z oi both sides of (^If), we get a relation between the Schwarz 
derivative and the r-function: 



nz) = Qz-' E 



z 



-k—n 



(9^ logr 



The latter can be used as an alternative definition of the r-function. 



(4^ 



7. Integrable structure of conformal maps. Eqs.( p^p^ ) allow one to say that the differ- 
ential 



dfl = Sdz + logw dtQ + '^{H^dtk — H^dtk) 

k=l 



(49) 



generates the set of Hamiltonian equations for deformations of the curve due to variation 
of tk. 

\S{z) = dA{z), di^S{z) = -d,Hk{z), (50) 

where we set Hq{z) = \ogw{z). The equations are consistent due to commutativity of 
the flows: 

Equations ( ^OD are more transparent being written in terms of canonical variables. The 
differential dfl suggests that the pairs logw, to and z{w), S{z{w)) are canonical and 
establishes the symplectic structure for conformal maps. Indeed, treating w as an inde- 
pendent variable, one rewrites eq. (B^ as 



{z{w), S{z{w))} 



(52) 
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where the Poisson bracket {, } is with respect to logw and the area to is defined as 



where the derivatives with respect to to are taken at fixed tk and w. 
The other fiows read 

^"^""^ -{H„z{w)} (54) 



dtk 
dS{z{w)) 



{Hu.S{z{w)}, (55) 



and similarly for the fiows with respect to t^- Now the Hamiltonian functions Hj. and 
H}, are degree k polynomials of w and respectively. 

The consistency conditions ( pO] ) now take the form of the zero- curvature conditions: 

dt^H,-dt^H, + {H,,H,} = Q, (56) 

dt^H, + d-t^H^ + H,} = 0. (57) 

The infinite set of the Poisson-commutating fiows form a Whitham integrable hierarchy 
|T0| . Eqs. (^55) are the Lax-Sato equations for the hierarchy. They generate an infinite 



set of differential equations for the coefficients (potentials) Uj of the inverse conformal 
map (0). The first equation of the hierarchy is 

dlt^^ = dt, exp{dtj), dtj = log r^. (58) 

The integrable hierarchy describing conformal maps is also known in the soliton literature 
as the dispersionless Toda lattice hierarchy, or SDiff(2) Toda hierarchy |Tl[] (see the 
next section)^ The algebra Sdiff(2) of area-preserving diffeomorphisms is the symmetry 
algebra of this hierarchy []1T|. Eqs. (^^-|57|) describe infinitesimal deformations of the 
curve such that the area to is kept fixed. 

The integrable hierarchy possesses many solutions. The particular solution relevant 
to conformal maps is selected by the subsidiary condition (|52D . This condition, known as 
dispersionless string equation, has already appeared in the study of the c = 1 topological 
gravity |]ll], |13|, |14| and in the large N limit of a model of normal random matrices [|15 
The latter is discussed in Sec. 9. 



8. Toda lattice hierarchy and its dispersionless limit. Below we review the two dimen- 
sional Toda lattice hierarchy and show that its dispersionless limit gives the equations 
describing conformal maps ( p5| , p6| , pl| , p5D . 

The 2D Toda hierarchy is defined by two Lax operators 

g OO g 

L = r{to) e^^—o + Mto) e"'^^ (59) 

k=0 

^ A relation between conformal maps of slit domains and special solutions to equations of hydrody- 
namic type (Benney equations) was first observed by Gibbons and Tsarev 
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L = e''—o r{to) + e'^'—oukito) (60) 

fc=0 

acting in the space of functions of to where the coefficients Uj and Uj are functions of to 
and also of two independent sets of parameters ("times") t^ and t^. Note that Uk and 
Uk as well as tk and t^ in (|59| , |60|) are not necessarily complex conjugate to each other, 
although we choose them to be so. 

The dependence of the coefficient Uk and Uk on t^ and t^ are given by the Lax-Sato 
equations: 

dL 

h— = [Hk,L] (61) 

Otk 

dL 

h-^ = [L,Hk] (62) 

and similar equations for L. The flows are generated by 

1 

2' 

where the symbol {L'')^ means positive (negative) parts of the series in the shift operator 
e ^'0 . The first equation of the hierarchy is the Toda lattice equation 

where r"^ = e</'(*o+?i)-0(io)_ 

The spectrum of the Lax operator is determined by the linear problem L\Ef = z"^. The 
wave function \E' is expressed through the r-function Tn of the dispersionfull hierarchy 
(^,|6^) by the following formula: 

^{z; to, ti, t2, ■ ■ ■) = rjT^to, ti, t2, . . .)z'°/''e^ E,>o'kz'' /^E,>o ^4ra(to, ti, ta, . . .) (66) 

Among many solutions of the hierarchy, one is of particular interest. It is selected by 



Hk = iL')^ + -{L% (63) 



Hk = m_ + -{L% (64) 



the string equation 16 



[L, L] = h (67) 

This solution is known to describe the normal matrix model at finite size of matrices [|15| . 

The dispersionless limit of the Toda hierarchy is a formal semi-classical limit h ^ 0. 

h-S- 

To proceed we notice that the shift operator W = e ^'0 obeys the commutation relation 
[W, to] = hW. In the semiclassical hmit it is supposed to be replaced by the canonical 
variable w with the Poisson bracket {logw,to} = 1. The Lax operator then becomes a 
c- valued function which is identified with the inverse conformal map z{w) (|2^) . Similarly, 
L is identified with S{z{w)). In their turn, the Lax-Sato equations ([61|J6^) are identified 



with eqs.(0,55) for the conformal map. In the same fashion the dispersionless limit 
of the string equation ( p7| ) is identified with eq. (|5^) . The semiclassical limits of the 
wave function and the r-function give the generating function Q and the dispersionless 
r-function: \E' — e^^^, r^ — >■ e*-'"^'^-'/^ . Similarly, eq. (|43|) is a semiclassical limit of the 
Hirota equation for the r-function of the 2D Toda hierarchy. 
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9. The T-function of the conformal map as large N matrix integral. The integrable 
structure of conformal maps is identical to the one observed in a class of random matrix 
models related to noncritical string theories. Moreover, there exists a random matrix 
model whose large limit reproduces exactly the r-function for analytic curves. 

Consider the partition function of the ensemble of normal random N x N matrices 
pffi, with the potential (|13D: 

Tn[t, t\= I ciMdMte-^^^(^'^') (68) 

A matrix is called normal if it commutes with its Hermitian conjugated [M, M^] = 0. 
Passing to the eigenvalues diag(2;i, . . . , z^) of the matrix M, one obtains the measure of 
the integral in a factorized form dMdM'^ ~ YiiLidzidzi Y[k<j{^k — Zj){zk — zj). Then 
the partition function is represents a two-dimensional Coulomb gas in the potential (|13D 

I- N 

rn[t,t\ = y n d^>^dzk e-^^^^^'^-^ j| g2iog|..-.,| (gg) 

fc=l i<j 

To proceed to the large limit one introduces a parameter to = 'hN and expresses the 
integrand in terms of density of eigenvalues as ^{py} where V^} is given by eq. 



(p!^). Then, the large N {h ^ 0) limit yields to the variational principle of Sec. 3. In the 
large A^ limit the eigenvalues of the matrix homogeneously fill the domain bound by 
the curve, characterized by the harmonic moments and the area to and leads to the 
T-function defined by eq. ([T6|) . Other objects introduced in Sees. 3-7 can also be identified 
with expectation values of the matrix model. In particular the moments Vk (eq.(^)) are 

Vk = n (ttm'') 

and Q' - ^vq = h (Tr log(^ - M)). 

In order to identify the Lax operator, we follow ||18|, [1^, |T^. Introduce the basis of 
orthogonal polynomials Pn{z) = hnZ"" + . . . (^ > 0), by the orthonormality relations 

(m|n) = [ dh Rjz) e-^^(^'^") P„(^) = (70) 



The polynomials are uniquely defined by the potential V up to phase factors. It is easy 
to see that the r-function is given by the product of the coefficients N\hnhn-i . . .ho of 
the highest powers of the polynomials Pn{z) = hnZ^ + ■ ■ ■. Then Lax operators L and 
L appear as the operators (jn\z\n) and {m\z\n). Since zPn{z) can be expressed through 
polynomials of the degree not grater than n, one may represent {m\z\n) and {m\z\n) in 

terms of shifts operators W = e ^'o in the form of ([59|J60|) , where r{to = hn) = hn/hn+i- 



Similar arguments allow one to identify the flows. Consider a variation of some 
operator {m\0\n) under a variation of t^. We have hdtf.{m\0\n) = {m[\Hk,0]\n), where 
Hk = Ak — A^k and {m\Ak\n) = {m\dt,.\n). Obviously Hk = —L''(W)+ negative powers of 
W. Choosing O to be Z (see (|60D) which consists on and positive powers of W, one 



concludes that if^ does not consists of negative powers of W. This brings us to eq . (|63D . 



^Earlier V. Kazakov pointed to us that the Lax equations (^j62) are generated by the Hermitian 
2-niatrix model JTzt with complex conjugated potentials. The latter and the normal matrix model have 
an identical 1/A^-expansion. 
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Finally, the operator D = {mlhdzln) is equal to 

D = L-J2ktkL''-\ (71) 

fc>i 

The Heisenberg relation [D,L] = h prompts the string equation (|57|). 

The matrix model also offers an effective method to derive eqs. ([37| - ^3D (see e.g. ||17|| ). 

10. T -function and spectral properties of the Dirichlet problem. This subject is under 
current study. 
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